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ABSTRACT 

A precision method for determining the value oi (3 = Q^/b, where b is the 
galaxy bias parameter, is presented. In contrast to other existing techniques 
that focus on estimating this quantity by measuring distortions in the redshift 
space galaxy-galaxy correlation function or power spectrum, this method removes 
the distortions by reconstructing the real space density field and determining 
the value of (3 that results in a symmetric signal. To remove the distortions, 
the method modifies the amplitudes of a Fourier plane-wave expansion of the 
survey data parameterized by /5. This technique is not dependent on the small- 
angle/plane-parallel approximation and can make full use of large redshift survey 
data. It has been tested using simulations with four different cosmologies and 
returns the value of (3 to ±0.031, over a factor of two improvement over existing 
techniques. 

Subject headings: cosmological parameters — galaxies: distances and redshifts — 
galaxies: statistics — large-scale structure of the universe — methods:data analysis 



1. Introduction 

One of the principal goals of the current large redshift surveys, such as the Two-Degree 
Field (2dF) and Sloan Digital Sky Surveys (SDSS), is an improved estimation of the total 
mass density parameter Q^n- A highly accurate measurement of Q^, along with measure- 
ments of the acoustic peaks in the CMB and the deceleration parameter from distant 
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supernovae, should enable cosmologists to pin down the global geometry of the Universe and 
infer the value of the vacuum energy density 

As is well-known, the redshift of a galaxy represents the sum of the galaxy's redshift 
distance, which depends upon cosmology, and its radial peculiar velocity. Therefore, raw 
redshift distance measurements are contaminated with the galaxies' radial peculiar veloci- 
ties. Consequently, galaxy redshift survey statistical measures such as the two-point galaxy 
correlation function contain distortions due to the existence of galaxy peculiar velocities. 

On large-scales in the linear regime, the galaxy peculiar velocities can be identified 
with fractional perturbations in the Hubble ratio AH/H due to fluctuations in the matter 
density field Ap/p. with AH/H oc Q^Ap/p} The peculiar velocities and resultant redshift 
distortions are statistical and by utilizing measures of the galaxy mass field, such as the 
galaxy-galaxy correlation function, together with the relation above, the distortions can be 
modeled and exploited as a way to measure ^2^- 

Since Qm is estimated using distortions in the galaxy correlation function or power 
spectrum and it is believed that galaxy formation is biased with respect to the underlying 
total mass field, what is actually measured is the parameter j3 = Q^;l^/b, where b is the 
bias in galaxy clustering with respect to the underlying mass field. Although somewhat 
problematic, this holds for any method that relies on the clustering of luminous matter to 
estimate a dynamical measure of the total mass field. 



1.1. Kaiser's Method 

The utility of using distortions in the galaxy redshift space correlation function to mea- 
sure P in the linear regime was developed in detail by Kaiser (1987), although its use was 
anticipated by others (see Sargent & Turner 1977; Peebles 1980 §76). In this seminal work. 
Kaiser's unique contribution was in showing that the anisotropics generated by these radial 
peculiar velocities could be modeled as redshift space density enhancements in the underlying 
mass field. 

Following Kaiser's lead, let us take a plane wave density fiuctuation 

A^(r) = AfeCOs(k-r + ^) (1) 

in real space. This real space plane wave fiuctuation maps to redshift space in the following 
simple way. 



^Assuming Ctx = 0. Otherwise O^*^ 0,%^ + %(1 -|- %^) (see Lahav et al. 1991) 
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A,(r)=A,(r)[l + 



1^ cos(k • f) ] 



(2) 



or 




line-of-sight f, is the redshift space density, and (3 — Jl^^/fe- In the hmit of small angular 
separations, the distortions in the redshift space power spectrum Ps(k) are given by 



where Pr(k) is the real space power spectrum. This limit is also known as the plane-parallel 
or distant observer approximation. Since the distortions contain powers of //^ and //^, they 
appear as a combination of quadrupole and hexadecapole distortions in Ps(k). As Kaiser 
went on to note, it seemed more practical to find a way to determine [3 by using this result 
to calculate distortions in the redshift space correlation function, the correlation function 
being the Fourier transform of the power spectrum. 



Traditionally, the approach to measuring distortions in the correlation function has 
been to decompose it in two-dimensions, where the two axes correspond to the directions 
parallel (tt) and perpendicular (rp) to the line-of-sight to a pair of galaxies. The line-of-sight 
direction has been taken to be the direction of the center-of-mass (Peebles 1980) or the half- 
angle (Landy, Szalay, & Broadhurst 1998; Hamilton 1992) between the pair. The resulting 
correlation function is denoted by tt). This decomposition for the correlation function 
for the simulation data is shown in Figure 1 on the following page. 

However, measuring j3 from distortions in the correlation function has proven difficult for 
several reasons. Firstly, as evident in Figure 1, the distortions at small separations arc due to 
galaxy pairwise velocities in the non-linear regime. These distortions are of an opposite sense 
than those due to large-scale hnear infall. Whereas the latter cause a compression in the 
correlation function along the line-of-sight, the non-linear distortions create an expansion. 
Therefore, they are competing effects and errors in their measurement may be positively 
correlated when both are estimated simultaneously. Secondly, each of these distortions affect 
the correlation function in terms of convolutions. Convolutions are troublesome since it 
is the unknown underlying function convolved with the kernel that is actually estimated. 
And, since the distortions enter in the radial direction only, calculations of the expected 



P,(k)=P,(k)(l + /?/.2)2 



(4) 



2. Difficulties in Measuring the Distortions 
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Fig. 1. — The real and redshift space correlation functions from one of the simulations. The 

expansion in the redshift space function at small separations along the n-axis is due to the 
small-scale non-linear pairwise velocity dispersion. The compression at larger separations is 
due to large-scale infall, which is a function of j3. 



correlation function given a particular power spectrum are extremely complicated and involve 
an integration over the spectrum itself, whose value may be only weakly constrained. For 
example, Kaiser's original paper derives the correlation function for the case of linear infall 
for just the line-of-sight direction as: 
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with 6=1. And thirdly, unless measurements of the correlation function are restricted 
to small angular separations for the galaxy pairs (the distant observer approximation), the 
expected distortions in the correlation function become even more complex. The complete 
problem in the case of wide-angles was recently solved in Szalay, Matsubara, & Landy (1998) 
and computer algebra was needed to keep track of the numerous terms. 

Attempting to measure the distortions using the power spectrum rather than the corre- 
lation function has also proven problematic. The basic intractability of this problem mainly 
results from the fact that radial distortions are not translationally invariant. The simple 
form for the distortions Ps(k) — Pr(is){l + /3/i^)^ given by Kaiser is only valid in the limit of 
small angles. Even with today's redshift surveys, this severely limits the amount of data that 
can be used in an analysis. A full decomposition of the distortions in the power spectrum 
in the case of wide-angles again results in very complex equations and requires a model for 
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the underlying power spectrum since the distortions cross-correlate Fourier modes on very 
different scales. 

Many researchers have been working over the last two decades on developing clever and 
robust methods to model the distortions and these methods have been based principally 
upon one of three approaches. The first involves measuring the ratio of the angle-averaged 

redshift space to real space power spectrum or correlation function (see Fry & Gaztaiiaga 
1994; Peacock & Dodds 1994; Baugh 1996; Loveday et al. 1996; Tadros & Efstathiou 1996; 
Peacock 1997). The second method calculates the ratio of the quadrupole to monopole 
moments of the redshift space power spectrum (see Hamilton 1993a,1995, 1997a; Bromley 
1994; Fisher et al. 1994; Cole, Fisher, & Weinberg 1994,1995; Lin 1995; Fisher & Nusser 
1996; Taylor & Hamilton 1996; Bromley, Warren, & Zurek 1997). The third utihzes a 
maximum likelihood approach in which the amplitudes of individual modes, /3, and the 
power spectrum are taken as parameters of the model (see Fisher, Scharf, & Lahav 1994; 
Heavens & Taylor 1995; Ballingcr, Heavens, & Taylor 1995). An excellent review of these 
methods along with a thorough background of the linear redshift distortion problem is given 
in Hamilton (1997a). 

Other work dealing primarily with the development of new methods, theoretical analysis, 
and constraints on surveys can be found in Lahav et al. (1991); Suto, & Suginohara (1991); 
Hamilton (1992,1993b,1997); Gramann, Cen, & Bachall (1993); Gramann, Gen, & Gott 
(1994); Hamilton & Gulhane (1996); Matsubara, & Suto (1996); Zaroubi, & Hoffman (1996); 
de Laix k Starkman (1998); Szalay Matsubara, & Landy (1998); Hatton & Gole (1998,1999); 
and Nakamura, Matsubara, & Suto (1998). 

In more recent work, Tadros et al. (1999) present a spherical harmonic analysis of the 
linear distortions in the PSGz Galaxy Gatalog and include a detailed synopsis of many of 
the linear distortion measurements listed above. They measure j3 = 0.58 ±0.26. In Matsub- 
ara, Szalay, & Landy (2000), the first application of the Karhunen-Love (K-L) eigenmode 
compression technique to this problem is presented using the Las Campanas Redshift Survey 
data with an estimation oi P — 0.30 ± 0.39. Hamilton, Tegmark, & Padmanabhan (2000) 
analyze the PSGz catalog also using a K-L expansion technique and determine P — 0.4llo;i2- 
A modified approach incorporating eigenvectors is described in Taylor et al. (2001) who find 
p = 0.39lo:}2- Ratcliffe et al. (1998) using the Durham/UKST Galaxy Redshift Survey 
estimate P ~ 0.5. Peacock et al. (2001), using over 141,000 galaxies from the Two-Degree 
Field Survey report a value oi P = 0.43 ± 0.07 using a two parameter model of the distor- 
tions in the redshift space correlation function. Tegmark, Hamilton, & Xu (2000) apply a 
pseudo K-L eigenmode analysis to the Two-Degree Field public release data and estimate 
/3 = 0.49 ± 0.16. 
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3. Removing the Distortions 

All of the methods listed above depend upon a measurement of the distortion signal, 
subject to all its complexities. However, Kaiser's identification of the redshift distortions with 
a statistically equivalent distortion of the real space mass density presents another approach 
to the problem. In Tegmark & Bromley (1995), an analytical reconstruction of a real space 
density field from a redshift space density field, such as that derived from a redshift survey, 
was described. However, as with attempts to model the distortions, the procedure had a 
complicated functional form and involved a relatively high dimensional integral. Because of 
this, it was suggested that such as inversion was best handled by parallel computation on 
a supercomputer. Setting computational complexities aside, the reconstruction of the real 
space density field presents a powerful new approach to this problem. 

As previously discussed, a direct fit of a model to either the correlation function or power 
spectrum in redshift space presents numerous difficulties due to the fundamental complexity 
of the effects of the distortions. Additionally, many of these methods include assumptions 
concerning either the underlying power spectrum or correlation function on large linear scales. 
However, rather than measuring the distortions themselves, another approach would be to 
invert the density enhancements and statistically re-create the real space density field. Since 
two-point estimators such as the power spectrum and correlation function are statistical 
anyway, one could invert the density field and simply search for the value of (5 that results 
in a symmetric function. In other words, when the field is properly inverted the anisotropies 
disappear. 

This method would work as follows: 

1) Take the data and perform the inversion for different values of the redshift distortions 

as parameterized by j3. 

2) For each of these inversions, calculate the power spectrum or correlation function 
that results from the inverted data. 

3) Look for the appropriate symmetry in the resulting function. 

Here the problem is reduced from the having to calculate and fit a very complex func- 
tion with an assumed power spectrum, to applying an inversion technique to the data itself 
and looking for a symmetric signal. 
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3.1. A Novel Method for Removing the Distortions 

Consider again Kaiser's analysis of the redshift space distortions. The power of his 
approach resuhed from considering a plane wave density fluctuation in real space. As shown 
above, this fluctuation maps to redshift space as 

A,(r)=A,(r)(l + /3/x2) (6) 

where r the line-of-sight in either real or redshift space. This relation is easily inverted to 
first order to give the real space density fluctuation, 



At first this development may seem academic in that redshift surveys do not appear to 
consist of plane wave density fluctuations. However, a redshift survey consists of a collection 
of galaxies, which for all intents and purposes can be considered point particles. As point 
particles, it is simple to construct a 5-function plane wave expansion of each individual point. 
Since each point has a well-defined position vector, it is easy to calculate the undistortion 
coefficient for each wavevector in the expansion for a given [5. Once this is done, 

the real space density can be reconstructed by re-summing this same set of plane waves with 
their modified amplitudes. 

Mathematically, the procedure is very straightforward. The real space mass density p(r) 
at any point r is given by 

Ngal n , ^ s 

^(r) = n 5] E (iT^j^-PfKkj • r - k, . rO]. (8) 

Here each galaxy i at position ri is expanded in a (5-function plane wave expansion over 
a set of n wavevectors kj, the amplitude of each wavevector is divided by the inversion 
coefficient, 1 -|- Pii^j where iiij = cos(kj • fi), the sum taken over all galaxies and the density 
reconstructed at r. Given the rotational but not translational invariance of these distortions, 
it is fortunate that the inversion of this problem can be expressed so simply in Cartesian 
coordinates using the Fourier transform. 

In practice this procedure is only slightly more complicated than that described above. 
The complications have to do with accounting for survey geometry, taking into account the 
non-linear small-scale pairwise distortions, choosing a finite set of wavevectors to expand 
over, and developing a sensitive method for determining when the (5 distortions have been 
removed. 



-8- 



4. Testing the Method 

To develop and test the method, numerical simulations were used (see Cole et al. 1999). 
Simulations are very convenient since the parameters of the 'data' are known and real space 
positions of the galaxies are also available. To make the procedure computationally efficient, 
the analysis was restricted to thin slices of simulated data 2° thick by 55° wide by 200 to 
500 h'^Mpc deep. Each slice contained approximately 5000 galaxies and the radial selection 
function was modeled after that of the SDSS. As an initial test, simulation data spanning 
four different cosmologies were used: two Cobe normalized models with f3 = 0.57 {flm = 0.5), 
/3 = 0.26 {Qm = 0.3), one tilted-flat model with f3 = 0.51 {flm = 1-0), and two independent 
rCDM fiat models with f3 = 0.54 (fi^ = 1.0). These are named 03, 05, E2, and E3SA, 
E3SB in Cole et al. (1999), respectively. 

Central to the method is to first expand the data in a Fourier expansion. Each point must 
be expanded and the new Fourier amplitudes calculated individually since the distortions 
depend upon cos(k ■ f). In computing the Fourier expansion, the data was binned with a 
resolution l/i^^Mpc using a grid of 1024 x 1024 in two dimensions and a FFT performed. By 
using very thin slices, the expansion could be carried out in two rather than three dimensions, 
limiting the computational effort. Further work will efficiently adapt this method to three 
dimensions. On this scale in two dimensions there are already 512 x 1024 amplitudes that 
must be calculated. 

After all amplitudes have been modified, their Fourier amplitudes are added back to- 
gether and an inverse FFT performed. This operation returns the undistorted density field. 
The field is then cut to correspond to the geometry of the original shce. This procedure was 
carried out for twenty slices and the Czi'^^p,'^) correlation function calculated for each one. 
All twenty correlation functions were averaged together to construct the mean. The mean 
correlation function was then windowed with a Hann window with a limit of 64/i~^Mpc and 
Fourier transformed. This returned the power spectrum of the Cz{'>^p,T^) decomposition. 

Although the focus of this method is to measure (3, as was shown in Figure 1, ^^(rp, tt) 
also contains small-scale pairwise velocity distortions in the non-linear regime. These dis- 
tortions have been shown to be well-characterized by an exponential distribution, which is 
equivalent to a Lorentzian distortion in the power spectrum (see Landy, Szalay & Broadhurst 
1998; Landy 2002). The full distortion function of the power spectrum is given by the ratio 

Psjk) ^ (1 + f3^'? . . 

Pr{k) l + \k^fj,^aUk) ^' 

where ai2/V^ is the decay width of the exponential. Using simulated data, it is simple to 
construct this ratio by dividing the power spectrum of the redshift space data by that of real 
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space data. The ratio of these power spectra is shown in Figure 2. The coordinate system 
is the Fourier analogue of the (r^, tt) basis and is denoted by {kr^, /Ctt)- 



Distortions in Simulation p = 0.57 IVIodel of p and Lorentzian Distortions 
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Fig. 2. — The left panel shows the distortions in the power spectrum for the simulated data. 
This was constructed by dividing the power spectrum in redshift space by the power spectrum 
in real space, using the real and redshift space functions from the simulated data. On the right 
is an analytical representation of the same function with j3 — 0.57 and (T12 = 360 km/sec. 
The agreement between the analytical model and the distortions is quite good. The 'bubble' 
in the center (long wavelengths) shows the linear infall distortions. The distortions across 
the graph are due to the small-scale pairwise velocity dispersion function (the Lorentzian). 

As shown in Figure 2, most of the signal from the /? distortions is found near the core. 
However, since the /? and small-scale distortions are competing effects, systematic errors will 
result if the small-scale distortions are not modeled correctly, especially at larger scales (in 
the core), and these errors will be positively correlated. A method to overcome the problem 
of modeling the small-scale distortions is described below. 

4.1. Removing the Distortions from Redshift Space Data 

To test the method, the procedure described above was used to remove the P distortions 
from the redshift space data. In practice, the distorted data was smoothed on 8 mpc scales 
and the correlation function windowed at 64 mpc to limit noise. Figure 3 shows the original 
smoothed and windowed power spectrum, and the power spectra for three values of (3 used 
to remove the distortions. 

Due to the residual non-linear distortions, the end result of this procedure is not expected 
to be a circularly symmetric function. However, the non-linear distortions effectively act only 
in the /c^ direction, that is along the line of sight. Therefore, even when the /3 distortions 
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Original Power Undistorted with Undistorted with Undistorted with 
13 = 0.57 13 = 0.35 |3 = 0.55 p = 0.75 
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Fig. 3. — The original power spectrum, and the power spectra for three sample values of (3 
used to remove the distortions for the j3 = 0.57 simulation. Although the case of j3 = 0.35 
appears to he the most symmetric result, the small scale non-linear distortions have not yet 
been taken into account. As will he shown in Figure 4, the vertically compressed signal of 
(5 — 0.55 is actually very close to the expected result. 

have been completely removed what is expected is a vertically compressed function about 
the krj, axis. 

What is needed at this stage is a method to account for the non-linear distortions. One 
of the fundamental challenges of any method is in trying to simultaneously fit the (5 and 
non-linear distortions. The main problem derives from the situation that these distortions 
are competing effects and their errors are positively correlated. 

To overcome this difficulty, a procedure was developed to estimate the non-linear distor- 
tions from the data itself, taking advantage of only the geometrical knowledge of how they 
effect the power spectrum. This also belies any problem which might arise from smoothing 
and windowing the data before transforming. 

In theory, neither the (5 nor the non-linear distortions effect the power along the k^^ 
axis, so the power along this axis can be utilized as the model underlying power spectrum. 
Additionally, if the (5 distortions have been removed in their entirety, the only remaining 
signal should be due to the residual non-linear distortions, which should be independent 

of krj,. Figure 4 shows the power spectrum of the data divided by the expected power 
as a function of |A:|, given by the signal along the k.f.^ axis. The power is almost perfectly 
symmetrical and independent of /c^p axis about the k,,.^ axis near the correct value of j3 = 0.57. 
Either under-correction using (3 = 0.35 or over-correction using j3 = 0.75, leave obvious 
signatures in the result. 
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Original Power 
p = 0.57 



Undistorted witli 
p = 0.35 



Undistorted witli 
p = 0.55 



Undistorted witli 
p = 0.75 




k/2n - Frequency Mpc"' 



Fig. 4. — The original power spectrum and that for three sample values of (3 used to remove 
the distortions divided by the model power spectrum as derived from the data (the value 
along the k^^ axis). If the /3 distortions have been correctly removed, this function should be 
symmetric about the axis. The cases for j3 = 0.35 and (3 = 0.75 clearly show under and 
over corrections. 



4.2. A Sensitive Fit 

The ratios shown in Figure 4 illustrate the utility of the technique. To actually deter- 
mine the best fit value of f3, however, the numerical difference between the intrinsic model 
derived from the data, that is the signal along the kr^ axis, and the undistorted spectra for 
different values of (3 was used. From this, a generalized 'chi-square' function was calculated 
to determine the optimal value of (3. Figure 5 shows this difference function for several values 
of (3 although over a much smaller range to illustrate the sensitivity of this method. 

|3 = 0.50 



(3 = 0.45 



13 = 0.55 |3 = 0.60 



o 



O 
0) 



CM 






k /27C - Frequency Mpc"^ 

Fig. 5. — The difference between four sample values of (3 and the underlying power derived 
from the signal along the k^-j, axis. Here a much smaller range in (3 was used in order to 
show the sensitivity of the method. The inversion of the differences is due to the different 
phase of under and over-corrections. 
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Figure 6 shows the 'chi-square' values for the appropriate P using four cosmologies and 
five simulations. This chi-squarc is simply the sum of the squares of the difference functions 
shown in Figure 5, normalized by dividing each by its minimum value for ease of comparison. 
The function was calculated inside a central ring of radius of /c/27r = 0.02, corresponding 
to wavelengths greater than or equal to 50 mpc. The sampling resolution of the difference 
function is AA;/27r = 0.002. 



Chi-Square Functions for Four Cosmologies - Five Simulations 



Open Model 03 
\ n = 0.30; |3 = 0.26 
\^ Measured p = 0.29 

\^ 30 


Open Model 05 . 50 
n = 0.50; p = 0.57 / 
, Measured p = 0.52 / 

\ / " 


Flat Model E2 

n = 1.0; p = 0.51 
Measured p = 0.48 


10 

8 

/ 4 


Flat Models E3S - A & B 
n = 1.0; p = 0.55 
Measured p^ = 0.55 ; 




\ ^ 10 



\ y 


2 







0.18 0.21 0.24 0.27 0.30 0.33 0.36 0.45 0.48 0.51 0.54 0.57 0.60 0.39 0.42 0.45 0.48 0.51 0.54 0.57 0.60 0.48 0.51 0.54 0.57 0.60 0.63 

Beta Inversion Value 



Fig. 6. — Chi-square functions for the method using CDM simulations with (3 = 0.26 (flm = 
0.3), i3 = 0.57 (Qm = 0.5), /S = 0.51 (Q„^ = 1.0), and f3 = 0.54 (Q^ = 1.0). These include 
COBE, structure normalized, and mixed models. The uncertainty in (3 across these results 
is ±0.031. For more detail on the simulations see the text and Cole et al. (1999). 

The results include a COBE normalized open /3 = 0.26 (Q^ = 0.3) model and a COBE 
normalized open (5 — 0.57 {Vim — 0.55) model, a mixed COBE and structure normalized 
fiat [3 = 0.51 {Qm = 1-0) model, and two independent structure normalized flat (3 = 0.54 
{Qrn = 1-0) models. Across these models, (3 has been measured to an uncertainty of ±0.031. 
The errors in the estimates arc not systematically high or low across the simulations. This 
compares favorably to the recent Two-Degree Field result of /3 = 0.43 ± 0.07 (see Peacock 
et al. 2001), where the uncertainties are over a factor of two higher. 



4.3. Difficulties in Applying the Method to Current Survey Data 

Although in theory this method is very straightforward, in practice the idiosyncracies of 
real survey data present other challenges. Fortunately, these difficulties will disappear once 
the surveys become more contiguous, as they result from the complexity of current survey 
geometries. 
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For the technique to be apphed successfully, it is necessary to be able to break up the 
survey into thin slices of fairly contiguous data, which do not include a large number of 
voids or holes. If there exist a large number of edges in the data due to a complex angular 
selection function, then the undistorted signal will leak out across the edges and into the 
geometrical voids in the survey when the data is inverted. This signal is then removed from 
the data when the geometric angular selection function cuts are made. This difficulty is not 
unexpected when considering that most of the distortion signal comes from scales above 30 
/i~^Mpc, and if the data is not contiguous on these scales, then the signal will be lost. To 
improve the method, experiments are underway to efficiently expand the method to three 
dimensions to circumvent this problem. 



5. Conclusion 

A method for determining the value oi (3 — /h based upon a Fourier inversion of 
distortions in the redshift space density field has been presented. In contrast to many other 
measures, this technique does not rely on a complex modeling of the expected distortions in 
the redshift space correlation function or power spectrum. Instead by inverting the redshift 
space density distortions as parameterized by it is possible to fit for a symmetric signal in 
the redshift space power spectrum. This technique is not dependent on the small angle/plane- 
parallel approximation and can make full use of large redshift survey data. The only present 
difficultly with this method is that it does depend on data that is fairly contiguous on 
large scales. Such data should become available from large redshift surveys within the 
next few years. It has been tested using simulations with four different cosmologies and 
returns the value of (5 to ±0.031, greater than a factor of two improvement over existing 
techniques. Presently, it has been tested in two-dimensions using thin slices of data to limit 
the computational burden, but can easily be expanded to three-dimensions. 
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